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We present new results for the properties of phases and phase transitions in spin-triplet ferromag- 
netic superconductors. The superconductivity of the mixed phase of coexistence of ferromagnetism 
and unconventional superconductivity is triggered by the presence of spontaneous magnetization. 
The mixed phase is stable but the other superconducting phases that usually exist in unconventional 
superconductors are either unstable or for particular values of the parameters of the theory some of 
them are metastable at relatively low temperatures in a quite narrow domain of the phase diagram. 
Phase transitions from the normal phase to the phase of coexistence is of first order while the phase 
transition from the ferromagnetic phase to the coexistence phase can be either of first or second 
order depending on the concrete substance. Cooper pair and crystal anisotropies determine a more 
precise outline of the phase diagram shape and reduce the degeneration of ground states of the 
system but they do not change drastically phase stability domains and thermodynamic properties 
of the respective phases. The results are discussed in view of application to metallic ferromagnets 
as UGe 2 , ZrZn 2 , URhGe. 

PACS numbers: 74.20.De, 74.20.Rp 



I. INTRODUCTION 

In 2000, experiments [l| at low temperatures (T~ IK) 
and high pressure (P ~ 1 GPa) demonstrated the exis- 
tence of spin triplet superconducting states in the metal- 
lic compound UGe2. The superconductivity is triggered 
by the spontaneous magnetization of the ferromagnetic 
phase that occurs at much higher temperatures. It coex- 
ists with the superconducting phase in the whole domain 
of its existence below T~ 1 K; see also experiments from 
Refs. and the discussion in Ref. Q. The same phe- 

nomenon of existence of superconductivity at low tem- 
peratures and high pressure in the domain of the (T, P) 
phase diagram where the ferromagnetic order is present 
was observed in other ferromagnetic metallic compounds 
(ZrZn 2 i and URhGe @) soon after the discovery Q of 
superconductivity in UGe2. 

In contrast to other superconducting materials, as 
ternary and Chevrel compounds, where the influence of 
magnetic order on superconductivity is also substantial 
(see, e.g., 0,011) El), in these ferromagnetic substances 
the phase transition temperature (T/) to the ferromag- 
netic state is much higher than the phase transition tem- 
perature (Tps) from ferromagnetic to a mixed state of co- 
existence of ferromagnetism and superconductivity. For 
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example, in UGe2, Tfs — 0.8 K while the critical temper- 
ature of the phase transition from paramagnetic to fer- 
romagnetic state in the same material is T) = 35K [10]. 
It can be assumed that in these substances the material 
parameter T s defined as the usual critical temperature 
of the second order phase transition from normal to uni- 
form (Meissner) supercondicting state in a zero external 
magnetic field is much lower than the phase transition 
temperature Tp$- The above mentioned experiments on 
the compounds UGe 2 , URhGe, and ZrZn 2 do not give 
any evidence for the existence of a standard normal-to- 
superconducting phase transition in a zero external mag- 
netic field. 

It seems that the superconductivity in the metallic 
compounds mentioned above always coexists with the fer- 
romagnetic order and is enhanced by it. In these systems, 
as claimed in Ref. the superconductivity probably 
arises from the same electrons that create the band mag- 
netism and can be most naturally understood rather as 
a triplet than spin-singlet pairing phenomenon. Metal- 
lic compounds UGc2, URhGe, and ZrZn 2 , are itinerant 
ferromagnets. An unconventional superconductivity is 
also suggested 11] as a possible outcome of recent exper- 
iments in Fe |l2j |. in which a superconducting phase has 
been discovered at temperatures below 2 K and pres- 
sures between 15 and 30 GPa. There both vortex and 
Meissner superconductivity phases 01 are found in the 
high-pressure crystal modification of Fe with a hexagonal 
close-packed lattice for which the strong ferromagnetism 
of the usual bec iron crystal probably disappears 0] . It 
can be hardly claimed that in hexagonal Fe the ferro- 
magnetism and superconductivity coexist but the clear 
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evidence for a superconductivity is also a remarkable 
achievement. 

The reasonable question whether these examples of 
superconductivity and coexistence of superconductivity 
and ferromagnetism are bulk or surface effects can be 
stated. The earlier experiments performed before 2004 do 
not answer this question. Recent experiments fl3| show 
that surface superconductivity appears in ZrZn 2 and its 
presence depends essentially on the way of preparation 
of the sample. But in our study it is important that bulk 
superconductivity can be considered well established in 
this substance. 

A phenomenological theory that explains the coexis- 
tence of ferromagnetism and unconventional spin-triplet 
superconductivity of Landau- Ginzburg type has been de- 
veloped recently in where possible low-order cou- 
plings between the superconducting and ferromagnetic 
order parameters are derived with the help of general 
symmetry group arguments. On this basis several im- 
portant features of the superconducting vortex state of 
unconventional ferromagnetic superconductors were es- 
tablished dm. 

In our paper we shall follow the approach from 
Refs. [3, [r| to investigate the conditions for the oc- 
currence of the Meissner phase and to demonstrate that 
the presence of ferromagnetic order enhances the p-wave 
superconductivity. We also establish the phase diagram 
of ferromagnetic superconductors in a zero external mag- 
netic field and show that the phase transition to the su- 
perconducting state can be either of first or second order 
depending on the particular substance. We confirm the 
predictions made in Refs. [bH ITU about the symmetry of 
the ordered phases. 

Our investigation is based on the mean-field approxi- 
mation |16| as well as on known results about the possible 
phases in non mag netic superconductors with triplet (p- 
wave) pairing F^flsL [TgL l2fl| . We extend our preceding 
results |2l], I22I |23| and show that taking into account the 
anisotropy of the spin-triplet Cooper pairs modifies but 
does not drastically change the thermodynamic proper- 
ties of the coexistence phase, especially in the tempera- 
ture domain above the superconducting critical temper- 
ature T s . The effect of crystal anisotropy is similar but 
we shall not make an overall thermodynamic analysis of 
this problem because we have to consider concrete sys- 
tems and crystal structures 0, |2(| for which there is no 
enough information from experiment to make conclusions 
about the parameters of the theory. Our results confirm 
the general concept that the anisotropy reduces the de- 
gree of ground state degeneration, and depending on the 
symmetry of the crystal, picks up a crystal direction for 
the ordering. 

There exists a formal similarity between the phase di- 
agram we obtain and the phase diag ram of certain im- 
proper ferroelectrics ||l El ElEl 0. We shall 
make use of the concept in the theory of improper fer- 
roelectrics, where the trigger of the primary order pa- 
rameter by a secondary order parameter (the electric po- 



larization) has been initially introduced and exploited; 
see Ref. [HEME! The mechanism of the M-triggered 
superconductivity in itinerant ferromagnets is formally 
identical to the mechanism of appearance of structural 
order triggered by the electric polarization in improper 
ferroelectrics. 

Our aim is to establish the uniform phases which are 
described by the GL free energy presented in Sec. II. Wc 
investigate a quite general GL model in a situation of a 
lack of concrete information about the values of the pa- 
rameters of this model for concrete compounds (UGe2, 
URhGc, ZrZn 2 ) where the ferromagnetic superconduc- 
tivity has been discovered. On the one hand the lack 
of information makes impossible a detailed comparison 
of the theory to available experimental data but on the 
other hand our results are not bound to one or more 
concrete substances and can be applied to any uncon- 
ventional ferromagnetic superconductor. In Sec. Ill the 
M-trigger effect will be described when only a linear cou- 
pling of the magnetization M to the superconducting 
order parameter %p is considered in a model of ferromag- 
netic superconductors where the spatial dependence of 
order parameters and all anisotropy effects are ignored. 
In Sec. IV we analyze the influence of quadratic coupling 
of magnetization to the superconducting order parameter 
on the thermodynamics of the ferromagnetic supercon- 
ductors. The application of our results to experimental 
(T, P) phase diagrams is discussed in Sec. IV. C. In Sec. 
V the anisotropy effects are outlined. In Sec. VI we 
summarize and discuss our findings. 



II. GINZBURG-LANDAU FREE ENERGY 

The general GL free energy functional, we shall use in 
our analysis, is 

F[^,M] = J d 3 xf(^,M), (1) 

where the free energy density f(ip, M) ( hereafter called 
"free energy") of a spin-triplet ferromagnetic supercon- 
ductor is a sum of five terms [lj, [15|, [17| , namely, 

f(i/>,M) = faW+&{M)+M<<l>,M)+%--B.M. (2) 

In Eq. (2) ip = {ipj;j — 1,2,3} is a three-dimensional 
complex vector describing the superconducting order and 
B = (H + 4-rrM) = V x A is the magnetic induction; H 
is the external magnetic field, A = {Aj;j = 1, 2, 3} is the 
magnetic vector potential. The last two terms on r.h.s. 
of Eq. (2) are related with the magnetic energy which 
includes both diamagnetic and paramagnetic effects in 
the superconductor; see, e.g., 0,H3- 

The term fs(ip) m Eq. (2) describes the superconduc- 
tivity for H = M = 0. It can be written in the form 
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Here 

/SW = f a ra d W+a s \M 2 + b fm 4 +Y l^| 2 + | E l^l 4 - 

(3) 



/sradW = ifi(A^)*(AVj) + ^2 [(AV>0*(A^) + (A^)*(A^)] + #3(AV>;)*(A^), (4) 

I 



where a summation over the indices «,j = 1,2,3 is as- 
sumed and the symbol 

of covariant differentiation is introduced. In Eq. (3), 
b s > and a s = a s (T — T s ), where a s is a positive 
material parameter and T s is the critical temperature of 
the standard second order phase transition which may 
occur at H = M = 0; H = \H\, and M = \M\. The 
parameters u s and v s describe the anisotropy of the spin- 
triplet Coo per pair and the crystal anisotropy, respec- 
tively, dill. Parameters Kj, (j = 1,2,3) in Eq. (4) 
are related with the effective mass tensor of anisotropic 
Cooper pairs [l7j . 

The superconducting part (3) of the free energy 
f(%j),M) is derived from symmetry group arguments and 
is independent of particular microscopic models; see, e.g., 
Refs. [13, El- According to classifications [lj, H3 the 
p-wave superconductivity in the cubic point group Oh 
can be realized through one-, two-, and three-dimensional 
representations of the order parameter. The expressions 
(3) and (5) incorporate all three possible cases. The co- 
efficients b s , u s , and v s in Eq. (3) are different for weak 
and strong spin-orbit couplings but in our investigation 
they are considered as undetermined material parameters 
which depend on the particular substance. 

The free energy of a standard isotropic ferromagnet is 
given by the term fp(M) in Eq. (2), 

MM) =c f J2 iViM.f + a f {T')M 2 + b -lM\ (6) 

where V, = d/dxj and 6/ > 0. The quantity <z/(Tx) = 
otf(T— Tj) is expressed by the material parameter a/ > 
and the temperature T'j which is different from the crit- 
ical temperature Tf of the ferromagnet and this point 
will be discussed below. We have already added a neg- 
ative term (— 2nM 2 ) to the total free energy f(tp,M) 
and that is obvious by setting H = in Eq. (2). The 
negative energy (— 2-itM 2 ) should be added to fp(M). 
In this way one obtains the total free energy f-p(M) of 
the ferromagnet in a zero external magnetic field that is 



given by a modification of Eq. (6) according to the rule 
f F (a f ) = [a f (T>) - a f (T f )} , (7) 
where a/ = a/(T — T/) and 



is the critical temperature of a standard ferromagnetic 
phase transition of second order. This scheme was used 
in studies of rare earth ternary compounds 0, H3, 0, 
l32l |. Alternatively [33|. one may use from the beginning 
the total ferromagnetic free energy /p(a/,Af) as given 
by Eqs. (6) - (8) but in this case the magnetic energy 
included in the last two terms on r.h.s. of Eq. (2) should 
be replaced with H 2 /8ir. Both approaches are equivalent. 
The term 

/i(^M)=i7 M.(^xf) + 5M 2 H 2 . (9) 

in Eq. (2) describes the interaction between the ferro- 
magnetic order parameter M and the superconducting 
order parameter if; [I4lll5| . The 70-term is the most sub- 
stantial for the description of experimentally found ferro- 
magnetic superconductors while the 5M 2 \ip\ 2 -term 
makes the model more realistic in the strong coupling 
limit as it gives the opportunity to enlarge the phase 
diagram including both positive and negative values of 
the parameter a s . In this way the domain of the stable 
ferromagnetic order is extended down to zero tempera- 
tures for a wide range of values of material parameters 
and the pressure P, a situation that corresponds to the 
experiments in ferromagnetic superconductors. 

In Eq. (9) the coupling constant 70 > can be rep- 
resented in the form 70 = AirJ, where J > is the 
ferromagnetic exchange parameter |15| . In general, the 
parameter S for ferromagnetic superconductors may take 
both positive and negative values. The values of the ma- 
terial parameters (T s , Tf, a s , ctf, b s , u s , v s , bf, Kj, 
70 and 5) depend on the choice of the concrete substance 
and on thermodynamic parameters as temperature T and 
pressure P. 

It is not easy to investigate straightforwardly the total 
free energy (2). In Ref. [l5j the authors used the crite- 
rion [34| for the stability of vortex state near the phase 



4 



transition line T C 2(H) (see also, Ref. J35|) and applied 
it with respect to the magnetization M. when H = for 
small values of \if)\ near the phase transition line T C 2(A4). 
We are interested in the uniform phases when the order 
parameters ip and M do not depend on the spatial vector 
x 6 V (V is the volume of the superconductor). There- 
fore, we present a detailed investigation of the coexis- 
tence of Meissner superconductivity and ferromagnetic 
order and, in particular, we show that the main prop- 
erties of the uniform phases can be described when the 
crystal anisotropy is ignored. We claim that some of the 
main features of the uniform phases in unconventional 
ferromagnetic superconductors can be reliably outlined 
even when the Cooper pair anisotropy is neglected. 

The magnetization M can be always assumed uniform 
outside a quite close vicinity of the magnetic phase transi- 
tion when the superconducting order parameter if) is also 
uniform, i.e., vortex phases are not present in the respec- 
tive temperature domain. These conditions are directly 
satisfied in type I superconductors but in type II su- 
perconductors the temperature should be sufficiently low 
and the external magnetic field should be zero. Never- 
theless, in type II superconductors these requirements for 
the appearance of uniform superconducting states may 
turn insufficient in materials having very high values of 
the spontaneous magnetization. In this case the uniform 
(Meissner) superconductivity and, hence, its coexistence 
with uniform ferromagnetic order may not occur even 
at zero temperature. Up to now type I unconventional 
ferromagnetic superconductors are not found experimen- 
tally. The predominant amount of experimental data for 
UGe2, URhGe, and ZrZn 2 do not give the possibility 
to conclude definitely either about the absence or the 
presence of uniform superconducting states at low and 
ultra-low temperatures but recently, an experimental ev- 
idence of uniform coexistence of superconductivity and 
ferromagnetism in UGe2 has been reported |36| . 

If real materials can be modelled by the general GL 
free energy (1) - (9), their ground state properties will 
be described by uniform states. The problem about the 
availability of such states in real materials at finite tem- 
peratures is quite subtle at the present stage of experi- 
mental research. We shall assume that uniform phases 
can exist in some unconventional ferromagnetic super- 
conductors, moreover these phases are solutions of the 
GL equations corresponding to the free energy (1) - (9). 
These arguments completely justify our study. 



In case of a strong easy axis type of magnetic 
anisotropy, as is in UGe2 p|, the overall complexity 
of mean-field analysis of the free energy f(if>, M) can 
be avoided by doing an Ising-like description: M = 
(0,0, M.), where M. = ±| JVfj is the magnetization along 
the z-axis. Because of the thermodynamic equivalence of 
up and down physical states (±Af) the analysis can be 
done only for M. > 0. But this approach can be also sup- 
ported without attracting crystal anisotropy arguments. 
When the symmetry of magnetic order is continuous, the 
symmetry of the total free energy f(if>, M) with respect 
to M comes into play and we can avoid the considera- 
tion of equivalent thermodynamic states that occur as a 
result of the respective symmetry breaking at the phase 
transition point but have no effect on thermodynamics of 
the system. In the isotropic system one may again choose 
the magnetization vector to point in the same direction 
as z-axis (\M\ = M z = M.) and this will not influence 
the generality of thermodynamic analysis. Here we pre- 
fer an alternative description for which the ferromagnetic 
state can occur as two thermodynamically equivalent up 
and down domains with magnetizations M. and (— M), 
respectively. 

We shall make the mean-field analysis of the uniform 
phases and the possible phase transitions between such 
phases in a zero external magnetic field (if = 0) when 
the crystal anisotropy is neglected (v s = 0) . The calcula- 
tions will be more easy to understand if we use notations 
that reduce the number of parameters in f(ip, M) by in- 
troducing 



b = (b s + u s + v s ). 



(10) 



Then we redefine the order parameters and all other pa- 
rameters in the following way: 

<Pj = b 1 / 4 ^ = 4> 3 e ie > , M = by 4 M , (11) 
a s a f u s v s 

r = Tb' t = Wf' W = T> 
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With the help of Eqs. (10) - (11) and using the uni- 
formity of if) and M we write the free energy density 
f(ifj, M) = F(if), M) /V, in the form 



f(if),M) = r0 2 + ^cf) 4 + 2j(f) 1 (f) 2 Msm(e 2 - 9^ + ll( f) 2 M 2 + tM 2 + ^M 4 (12) 



-2w [</> 2 </> 2 sin 2 (0 2 - 8 X ) + 0?0§sin 2 (0i - 3 ) + ^I;sin 2 (02 - 3 )] " v[fy 2 2 + 4>\4>l + 4>l^, 



In the above expression the order parameters if) and M are defined per unit volume. 
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The equilibrium phases are obtained from the equa- 
tions of state 



df(po) 



0. 



(13) 



where [i = {^ a } = (M, (pi, </> 3 , 9i,...,9 3 ) and fi de- 
notes an equilibrium phase. The stability matrix F of 
the phases /j,q is given by 



d 2 /(/io) 



(14) 



An alternative treatment can be done in terms of real 
(?/>j) and imaginary (V 1 ") parts of the complex numbers 
ijjj = ip'j + iipj . The calculation with moduli <pj and phase 
angles 9j of ipj is more simple but in cases of strongly de- 
generate phases some of the angles 9j remain unspecified. 
Then an alternative analysis with the help of the compo- 
nents tjj'j and tp'j should be done. 

The thermodynamic stability of the phases that are 
solutions of Eqs. (13) is checked with the help of the ma- 
trix (14). An additional stability analysis is done by the 
comparison of free energies of phases that satisfy (13) and 
render the stability matrix (14) positive in one and the 
same domain of parameters {r, t, 7, 71, w, v }. This step is 
important because the complicated form of the free en- 
ergy generates a great number of solutions of Eqs. (13) 
and we have to sift out the stable from metastable phases 
that correspond either to global or local minima of the 
free energy, respectively |l6j . 

Some solutions of Eqs. (13) have a marginal stability, 
i.e., their stability matrix (14) is neither positively nor 
negatively definite. This is often a result of the degener- 
ation of phases with broken continuous symmetry. If the 
reason for the lack of a clear positive definiteness of the 
stability matrix is precisely the mentioned degeneration 
of the ground state, one may reliably conclude that the 
respective phase is stable. If there is another reason, the 
analysis of the matrix (14) will be insufficient to deter- 
mine the respective stability property. These cases are 
quite rare and occur for particular values of the parame- 
ters {r,f,7, ...}. 



III. SIMPLE CASE OF M-TRIGGERED 
SUPERCONDUCTIVITY 

We shall consider the Walker-Samokhin model pl| 
when only the M(j>\ 4>2— coupling between the order pa- 
rameters i/j and M is taken into account (7 > 0, 71 = 0) 
and the anisotropics (w = v = 0) are ignored. The uni- 
form phases and the phase diagram in this case were in- 
vestigated in Refs. [21], |22, H3 ■ Here we summarize the 
main results in order to make a clear comparison with 
the new results presented in Sections IV and V. Our 
main aim is the description of a trigger effect which con- 
sists of the appearance of a "compelled superconductiv- 
ity" caused by the presence of ferromagnetic order (here, 



this is a standard uniform ferromagnetic order); see also 
Refs. pll I22I I23I ] where this effect has been already estab- 
lished and briefly discussed. As mentioned in the Intro- 
duction, a similar trigger effect is known in the physics 
of improper ferroelectrics. We shall set #3 = and use 
the notation 9 = A9 = (9 2 - di). 



A. Phases 

The possible (stable, metastable and unstable) phases 
are given in Table 1 together with the respective existence 
and stability conditions. The normal or disordered phase, 
denoted in Table 1 by N, always exists (for all tempera- 
tures T > 0) and is stable for t > 0, r > 0. The super- 
conducting phase denoted in Table 1 by SCI is unstable. 
The same is valid for the phase of coexistence of ferro- 
magnetism and superconductivity denoted in Table 1 by 
C02. The N-phase, the ferromagnetic phase (FM), the 
superconducting phases (SC1-3) and two of the phases 
of coexistence (COl-3) are generic phases because they 
appear also in the decoupled case (7 = 0). When the 
M(f>x 02-coupling is not present, the phases SC1-3 are 
identical and represented by the order parameter cj> with 
components tfij that participate on equal footing. The 
asterisk attached to the stability condition of the second 
superconductivity phase (SC2) indicates that our analy- 
sis is insufficient to determine whether this phase corre- 
sponds to a minimum of the free energy. It will be shown 
that the phase SC2, two other purely superconducting 
phases and the coexistence phase COl, have no chance 
to become stable for 7 7^ 0. This is so, because the phase 
of coexistence of superconductivity and ferromagnetism 
(FS in Table 1), that does not occur for 7 = 0, is stable 
and has a lower free energy in their domain of stability. 
A second domain (M < 0) of the FS phase exists and is 
denoted in Table 1 by FS*. Here we shall describe only 
the first domain FS. The domain FS* is considered in the 
same way. 

The cubic equation for magnetization of FS-phase (see 
Table 1) is shown in Fig. 1 for 7 = 1.2 and t = —0.2. For 
any 7 > and t, the stable FS thermodynamic states are 
given by r(M) < r m = r(M m ) for M > M m > 0, where 
M m corresponds to the maximum of the function r(M). 
The dependance of M m (t) and M (f) = (-t + j 2 /2) 1 / 2 = 
\f~ZM m {t) on t is drawn in Fig. 2 for 7 = 1.2. Functions 
r m (t) = 4M^(t)/ 7 for t < j 2 /2 (depicted by the line of 
circles in Fig. 3) and 



r e (t)= 7 |*| 1/2 , 



(15) 



for t < define the borderlines of stability and existence 
of FS. 



B. Phase diagram 

We have outlined the domain in the (i, r) plane where 
the FS phase exists and is a minimum of the free en- 




FIG. 1: A = 7r/2 as a function of M for 7 = 1.2, and t = 



—0.2. The parameters r, t, and 7 are given by Eq. (11). 



7 




FIG. 2: The magnetization M versus t for 7 = 1.2: the dashed line represents M , the solid line represents M eq , and the dotted 
line corresponds to M m . 
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FIG. 3: The phase diagram in the plane (t, r) with two tricritical points (A and B) and a triple point C; 7 = 1.2. The 
parameters r ~ [T — T a (P)] and t ~ [T — T/(P)] are defined by Eq. (11). The domains of existence and stability of the phases 
N, FM and FS are shown. The line of circles represents the function r m (t) given by Eq. (17). The dotted line represents the 
function r e (f) given by Eq. (15). On the left of point B, the same dotted curve corresponds to a FM-FS phase transition of 
second order. The equilibrium lines of N-FS and FM-FS phase transitions of first order are given by the solid lines AC and 
CB, respectively. 
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TABLE I: Phases and their existence and stability properties [9 = [Qi — 9i), k — 0, ±1, ...]. 



Phase 


order parameter 


existence conditions 


stability domain 


IN 


(pj = Ivl = U 


always 


t > U, r > U 


FM 


<j>j = 0, M 1 = -t 


t < 


r > 0, r > r e (t) 


SCI 


J. l\ /T n J.2 

<pi = M = U, = — r 


r < 


unstable 


SC2 


0^ = -r, = Trfc, M = 


r < 


(t > 0)* 


SC3 


0i = 02 = M = 0, 4>'i = -r 


r < 


r < 0, t > 


COl 


0i = 02 = 0, <f>i = -r, = -t 


r < 0, t < 


r < 0, t < 


C02 


0i = 0, 0^ = -r, 6 = 62 = irk, M z = -t 


r < 0, t < 


unstable 


FS 


24>{ = 20^ = 0^ = -r + jM, 03 = 

6> = 27r(fc - 1/4), 7 r = (7 2 - 2t)M - 2M 3 


7M > r 


3M' Z > (-t + 772) 
M > 


FS* 


20? = 20^ = 0^ = -0 + jM), 03 = 
= 2n(k + 1/4), 7r = (2t - -y 2 )M + 2M 3 


-7M > r 


3M' Z > (-t + ^/2) 
M < 



ergy. For r < the cubic equation for M (see Table 
1) and the existence and stability conditions are satis- 
fied for any M > provided t > 7 2 . For t < 7 2 the 
condition M > Mo have to be fulfilled, here the value 
M = (-t + 7 2 /2) 1/2 of M is obtained from r(M ) = 0. 
Thus for r = the N-phase is stable for £ > 7 2 /2, and 
FS is stable for t < 7 2 /2. For r > 0, the requirement for 
the stability of FS leads to the inequalities 



max I — , M„ 

,7 



< M < Mr 



0, 



(16) 



where M m = (Mo/y/S) and M should be the positive 
solution of the cubic equation of state from Table 1; 
M m > gives a maximum of the function r(M); see 
also Figs. 1 and 2. 

The further analysis defines the existence and stability 
domain of FS below the line AB denoted by circles (see 
Fig. 3). In Fig. 3 the curve of circles starts from the 
point A with coordinates ("f 2 /2, 0) and touches two other 
(solid and dotted) curves at the point B with coordinates 
(ie = — 7 2 /4, tb = 7 2 /2). Line of circles represents the 
function r(M m ) = r m (t) where 



4 (1 2 
3^7 V 2 



t 



3/2 



(17) 



Dotted line represents r e (t), defined by Eq. (15). The 
inequality r < r m (t) is a condition for the stability of 
FS, whereas the inequality r < r e (t) for (— t) > 7 2 /4 is 
a condition for the existence of FS as a solution of the 
respective equation of state. This existence condition for 
FS is obtained from 7M > r (see Table 1). 

In the region on the left of the point B in Fig. 3, the 
FS phase satisfies the existence condition 7M > r only 
below the dotted line. In the domain confined between 
the lines of circles and the dotted line on the left of the 
point B the stability condition for FS is satisfied but the 
existence condition is broken. The inequality r > r e (t) is 
the stability condition of FM for < (-i) < 7 2 /4. For 
(— t) > 7 2 /4 the FM phase is stable for all r > r e (t). 

In the region confined by the line of circles AB, the 
dotted line for < (— t) < 7 2 /4, and the t— axis, the 



phases N, FS and FM have an overlap of stability do- 
mains. The same is valid for FS, the SC phases and 
COl in the third quadrant of the plane (t, r). The com- 
parison of the respective free energies for r < shows 
that the stable phase is FS whereas the other phases are 
metastable within their domains of stability. 

The part of the i-axis given by r — and t > 7 2 /2 is a 
phase transition line of second order which describes the 
N-FS transition. The same transition for < t < 7 2 /2 
is represented by the solid line AC which is the equilib- 
rium transition line of a first order phase transition. The 
equilibrium transition curve is given by the function 



,(*) 



3 7 - ( 7 2 + 16t) 1/2 M eq (t). (18) 



Here 



M eq (t) - ^-j= [ 7 2 - 8t + 7 ( 7 2 + 16t) 



1/2 



1/2 



(19) 



is the equilibrium jump of the magnetization. The order 
of the N-FS transition changes at the tricritical point A. 

The domain above the solid line AC and below the line 
of circles for t > is the region of a possible overheating 
of FS. The domain of overcooling of the N-phase is con- 
fined by the solid line AC and the axes (t > 0, r > 0). 
At the triple point C with coordinates [0, r eq (0) — 7 2 /4] 
the phases N, FM, and FS coexist. For t < the straight 
line 



< q {t) = \ + \t\, 



t B <t < 0, 



(20) 



describes the extension of the equilibrium phase transi- 
tion line of the N-FS first order transition to negative 
values of t. For t < ts the equilibrium phase transition 
FM-FS is of second order and is given by the dotted line 
on the left of the point B which is the second tricriti- 
cal point in this phase diagram. Along the first order 
transition line r* q (t) given by Eq. H2U|) the equilibrium 
value of M is M eq =7/2, which implies an equilibrium 
order parameter jump at the FM-FS transition equal to 
(7/2 — x/Jtf). On the dotted line of the second order FM- 
FS transition the equilibrium value of M is equal to that 
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of the FM phase (M eq = y/\t\). The FM phase does not 
exist below T s and this is a shortcoming of the model 112(1 
with 71 = 0. 

The equilibrium FM-FS and N-FS phase transition 
lines in Fig. 3 can be expressed by the respective equi- 
librium phase transition temperatures T eq defined by the 
equations r e = r(T eq ), r eq = r(T eq ), r* eq = r(T eq ), and 
with the help of the relation M eq — M(T eq ). This limits 
the possible variations of parameters of the theory. For 
example, the critical temperature (T eq = T c ) of the FM- 
FS second order transition (7 2 /4 < —t) is obtained in the 
formT c = (T s + 4:nJM/a s ), or, using M = (-0//6/ ) 1/2 , 



T r = T s - 



T*{T f -T s ) 



1/2 



(21) 



Here Tf > T s , and T* = (47rJ) 2 a//a 2 6/ is a characteris- 
tic temperature of the model 1(12(1 with 71 = w = v = 0. 
A discussion of Eq. (21) is given in Sec. IV. C. 

The investigation of the conditions for the validity of 
Eq. l(2*T|l leads to the conclusion that the FM-FS contin- 
uous phase transition (at 7 2 < —t) will be possible only 
if the following condition is satisfied: 



T f -T s >=(<; + V?)T* 



(22) 



where <r = &/Q; 2 /46 s a 2 . Therefore, the second order FM- 
FS transition should disappear for a sufficiently large 7- 
coupling. Such a condition does not exist for the first 
order transitions FM-FS and N-FS. 

The inclusion of the gradient term (4) in the free 
energy J5J should lead to a depression of the equilib- 
rium transition temperature. As the magnetization in- 
creases with the decrease of the temperature, the vor- 
tex state should occur at temperatures which are lower 
than the equilibrium temperature T eq of the Meissner 
state. For example, the critical temperature (T c ) corre- 
sponding to the vortex phase of FS-type has been evalu- 
ated [l5j to be lower than the critical temperature 1(21(1 : 
(T c — T c ) = Air^sM/ccs, where fis = \e\h/2mc is the 
Bohr magneton. For J fis, we have T c w T c . 

For r > 0, namely, for temperatures T > T s the super- 
conductivity is triggered by the magnetic order through 
the 7-coupling. The superconducting phase for T > T s is 
entirely in the (t, r) domain of the ferromagnetic phase. 
Therefore, the uniform supeconducting phase can occur 
for T > T s only through a coexistence with the ferromag- 
netic order. 

In the next Sections we shall focus on the temperature 
range T > T s which seems to be of main practical inter- 
est. We shall not dwell on the superconductivity in the 
fourth quadrant [t > 0, r < 0) of the (t, r) diagram where 
pure superconducting phases can occur for systems with 
T s > Tf, but this is not the case for UGe2, URhGe and 
ZrZn 2 . Also we shall not discuss the possible metastable 
phases in the third quadrant (t < 0,r < 0) of the (t, r) 
diagram. 



C. Magnetic susceptibility 

We consider the longitudinal magnetic susceptibility 
Xi = (xv/V) P er un it volume j^. The external mag- 
netic field H = (0,0, H) with H = (df/dM) has the 
same direction as the magnetization M. We shall calcu- 
late the quantity \ = y/bfXi f° r the equilibrium thermo- 
dynamic states fio given by Eq. (13). Having in mind the 
relations (11) between M and Ai, and between i/j and ip 
we can write 



dM 



df_\ 
dM) 



(23) 



where the equilibrium magnetization Mo and equilibrium 
superconducting order parameter components tpoj should 
be taken for the respective equilibrium phase. See Ta- 
ble 1, where the suffix "0" of cf>, 8, and M is omitted; 
hereafter this suffix will be often omitted. The value of 
the equilibrium magnetization M in FS is the maximal 
nonnegative root of the cubic equation in M given in 
Table 1. 

From Eq. (23) we obtain the susceptibility \ of FS 
phase in the form 



X 



-7 2 + 2t + 6M 2 



(24) 



The susceptibility of the other phases has the usual ex- 
pression 



X- 1 =2t 



6M Z 



(25) 



Eq. (25) yields as results the paramagnetic susceptibility 
(\p = l/2t; t > 0) of the normal phase and the fer- 
romagnetic susceptibility (\f = 1/4 j tj ; t < 0) of FM. 
These susceptibilities can be compared with the suscep- 
tibility x of FS which cannot be calculated analytically in 
the whole domain of stability of FS. Therefore, we shall 
consider the close vicinity of the N-FS and FM-FS phase 
transition lines. 

Near the second order phase transition line on the left 
of the point B (t < ts), the magnetization has a smooth 
behavior and the magnetic susceptibility does not exhibit 
any singularities like jump or divergence. For t > 7 2 /2, 
the magnetization is given by M = (s_ + s+), where 



s± = 



' 4 



27 



1/2' 



1/3 



(26) 



When r = 0, it is obvious that also M — 0. For \^r\ <C 
(t - 7 2 /2) we have M w -jr/(2t - 7 2 ) < 2t. Therefore, 
in a close vicinity (r < 0) of r = along the second 
order phase transition line (r = 0, t > 7 2 /2) the magnetic 
susceptibility is well described by the paramagnetic law 
Xp = (1/2*). For r < and t -> 7 2 /2, we obtain M = 
— (7r/2) 1 / 3 which gives 



X- 1 =6 



j\r\ 



2/3 



(27) 
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On the phase transition line AC 
1 



M eq {t) 



and, hence, 



X 



2^2 



7 2 - 8t + 7 (7 2 + 16t) 



1/2 



1/2 



-4t- 



1-3 1 



16t 



1/2' 



(28) 



(29) 



At the tricritical point A this result gives \ = 0, 

and at the triple point C with coordinates (0, 7 2 /4) we 
have x(C) = (2/7 2 ). On the line BC we obtain M = 7 /2 



so 



X 



(30) 



At the tricritical point B with coordinates (— 7 2 /4, 7 2 /2) 
the result is x _1 (.B) = 0. 

To investigate the magnetic susceptibility tensor we 
shall consider arbitrary orientations of the vectors H and 
M. We denote the spatial directions (as, y, z) by (1, 2, 3). 

The components of the inverse magnetic susceptibility 
tensor 



Xi 



r 1 Vb~f = {Xu 1 }Vb~f 



(31) 



can be represented in the form 



<:) 



2(t + M 4 )Sn + AMiMj + il-Q^iv X ¥>*)<, (32) 



where M and tpj are taken at their equilibrium values: 
Mo, </?oj, 8oj. The last term in r.h.s. of Eq. (29) is equal 
to zero for all phases in Table 1 except for FS and FS*. 
When the second term in Eq. (30) is equal to zero we 
obtain the known result of the susceptibility tensor for 
second order phase transitions; see, e.g., |16|. 

In FS phase <f)j depend on Mj and we can choose again 
M = (0, 0, M) and use the results from Table 1 for the 
equilibrium values of 9 and M. Then the components 
Xij 1 corresponding to FS are 



Xil = 2(t + M 2 )5ij + AMiMj - 7 2 ^ 
Thus we have X&j = 0; 



(33) 



(34) 



and X33 coincides with the inverse longitudinal suscepti- 
bility x _1 as given by Eq. (24). 



D. Entropy and specific heat 

The entropy S(T) = (S/V) = -Vd(f/dT) and the 
specific heat C(T) = (C/V) = T(dS/dT) per unit vol- 
ume V are calculated in a standard way 161. We are 
interested in the jumps of these quantities on the N-FM, 



FM-FS, and N-FS transition lines. The behavior of S(T) 
and C(T) near the N-FM phase transition and near the 
FM-FS phase transition line of second order on the left of 
the point B (Fig. 3) is known from the standard theory 
of critical phenomena and for this reason we focus our 
attention on the first order phase transitions FS-FM and 
FS-N for t > — 7 2 /4, i.e., on the right of the point B in 
Fig. 3. 

We make use of the equations for the order parameters 
tp and M from Table 1 and apply the standard procedure 
for the calculation of S: 



S(T) = 



=M 2 . 



(35) 



The next step is to calculate the entropies Sfs{T) and 
S f m of the ordered phases FS and FM. We shall stick to 
the usual convention Fm = V/jv = for the free energy 
of the N-phase , so we must set 5jv(T) = 0. 

Near the second order phase transition line (r = 0, 
t > 7 2 /2), Sfs(T) is a smooth function of T and has no 
jump but the specific heat Cfs has a jump at T = T s , 
i.e. for r = 0. This jump is given by 



AC FS (T S ) 



1 



1 



1 - 2t(T s )/ 7 2 



(36) 



The jump ACfs(T s ) is higher than the usual jump 
AC(T C ) = T c a 2 s /b s known from the Landau theory of 
standard second order phase transitions |16| . 

The entropy jump AS A c(T) = S F s(T) on the line AC 



AS AC (T) 



(37) 



-M 



eg 



4x/67 



16* 

„,2 



1/2- 



where M eq is given by Eq. (19). From Eqs. (19) and 
(37), we have AS(t = 7 2 /2) = 0, i.e., AS(T) becomes 
equal to zero at the tricritical point A. We find also from 
Eqs. (19) and (37) that at the triple point C the entropy 
jump is 




(38) 




On the line JSC the entropy jump is defined by 
AS BC {T) = [S FS {T) - S FM {T)}. We obtain 



(39) 



At the tricritical point B this jump is equal to zero as 
should be. The calculation of the specific heat jump on 
the first order phase transition lines AC and BC is re- 
dundant for two reasons. Firstly, the jump of the specific 
heat at a first order phase transition differs from the en- 
tropy by a factor of order of unity. Secondly, in caloric 
experiments where the relevant quantity is the latent heat 
Q = TAS(T), the specific heat jump can hardly be dis- 
tinguished. 
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E. Note about a simplified theory 

The analysis in this Section can be done following an 
approximate scheme known from the theory of improper 
ferroelectrics; see, e.g., Ref. [2^. In this approxima- 
tion the order parameter M is considered small enough 
which makes possible to ignore M 4 -term in the free en- 
ergy. Then one easily obtains from the data for FS 
presented in Table 1 or by a direct calculation of the 
respective reduced free energy that the order parame- 
ters 4> an d M of FS-phase are described by the simple 
equalities r = (7M - 2 ) and M = (7/2i)0 2 . For fer- 
roelectrics working with oversimplified free energy give s 
a substantial departure of theory from experiment |29| . 
For ferromagnetic superconductors the domain of relia- 
bility of this approximation could be the close vicinity of 
the ferromagnetic phase transition, i.e. for temperatures 
near the critical temperature T/ . This discussion can be 
worthwhile if only the primary order parameter also ex- 
ists in the same narrow temperature domain (0 > 0). 
Therefore, the application of the simplified scheme can 
be useful in systems, where T s > Tt. 

For T s < Tf , the analysis can be simplified if we sup- 
pose a relatively small value of the modulus of the 
superconducting order parameter. This approximation 
should be valid in some narrow temperature domain near 
the line of second order phase transition from FM to FS. 



IV. EFFECT OF SYMMETRY CONSERVING 
COUPLING 

Here we shall include in our consideration both lin- 
ear and quadratic couplings of magnetization to the su- 
perconducting order parameter which means that both 
parameters 7 and 71 in free energy (12) are different 
from zero. In this way we shall investigate the effect 
of the symmetry conserving 71-term in the free energy 
on the thermodynamics of the system. When 7 is equal 
to zero but 71 7^ the analysis is easy and the results 
are known from the theory of bicritical and tetracritical 
points 0, 127], H3, . For the problem of coexistence of 
conventional superconductivity and ferromagnetic order 
the analysis (7 = 0, 71 7^ 0) was made in Ref. 

At this stage we shall not take into account any 
anisotropy effects because we do not want to obscure 
the influence of quadratic interaction by considering too 
many parameters. For 7, 71 7^ and w = 0, v = the 
results again can be presented in an analytical form, only 
a small part of phase diagram should be calculated nu- 
merically. 



A. Phases 

The calculations show that for temperatures T > T s , 
i.e., for r > 0, we have again three stable phases. Two 
of them are quite simple: the normal (N-) phase with 



existence and stability domains shown in Table 1, and the 
FM phase with the existence condition t < as shown in 
Table 1, and a stability domain defined by the inequality 



A 1 ^ < r. Here 



7i* + 7 



(40) 



and one can compare it with the respective expres- 
sion (15) for 71 = 0. In this paragraph we shall retain 
the same notations as in Sec. Ill, but with a superscript 
(1) in order to distinguish them from the case 71 = The 
third stable phase for r > is a more complex variant 
of the mixed phase FS and its domain FS*, discussed in 
Sec. III. The symmetry of the FS phase coincides with 
that found in fl5| . 

We have to mention that for r < there are five pure 
superconducting (M = 0, > 0) phases. Two of them, 
(0i > 0,02 = 03 = 0) and (0i = O,0 2 > 0,03 > 0) 
are unstable. Two other phases, (0i > 0, 02 > 0, 03 = 



0,02 



6»i + 7rfc) and (0i > 0, 2 > 0, 3 > 0, 1 



9\ + 7rfc, $3 - arbitrary; k — 0,±1,...) show a marginal 
stability for t > 71 r. 

Only one of the five pure superconducting phases, the 
phase SC3, given in Table 1, is stable. In case of 71 7^ 
the values of 4>j and the existence domain of SC3 are the 
same as shown in Table 1 for 71 = but the stabil- 
ity domain is different and is given by t > "fir. When 
the anisotropy effects are taken into account the phases 
exhibiting marginal stability within the present approx- 
imation may become stable. Besides, three other mixed 
phases (M 7^ 0, > 0) exist for r < but one of them 
is metastable (for 7 2 > l,t < 71 r, and r < 71 1) and 
the other two are absolutely unstable. Here the ther- 
modynamic behavior for r < is much more abundant 
in phases than for improper ferroelectrics with two com- 
ponent primary order parameter |27| . However, at this 
stage of experimental needs about the properties of un- 
conventional ferromagnetic superconductors the investi- 
gation of the phases for temperatures T < T s is not of 
primary interest and for this reason we shall focus our 
attention on the temperature domain r > 0. 

The FS phase for 71 7^ is described by the following 
equations: 



= fa = -j= , 03 = 0, 



2 = (±jM -r -7iM 2 ), 



(41) 



(42) 



(l- 7l 2 )M 3 ±^77iAf 2 + (t-^L 



and 



(0 2 



7l rjM±^ = 0, 
(43) 

27rfc, (44) 
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(k = 0, ±1, ...)■ The upper sign in Eqs. (42) - (44) cor- 
responds to the FS domain where sin(0 2 — #i) = — 1 
and the lower sign corresponds to the FS* domain with 
sin(6> 2 — 6*i ) = 1. This is a generalization of the two- 
domain FS phase discussed in Sec. III. The analysis of 
the stability matrix (14) for these phase domains shows 
that FS is stable for M > and FS* is stable for M < 0, 
just like our result in Sec. III. As these domains belong 
to the same phase, namely, have the same free energy 
and are thermodynamically equivalent, we shall consider 
one of them, for example, FS. 



r (i) 

eq 



(*) 



(48) 



M. 



eg 



( 7l M eg - 7 /2) 



(1 - j 2 )M 2 eq + | 7 7iM e9 + (t~\) 



and to substitute the above expression in the existence 
and stability conditions of FS-phase. It is obvious that 
there is a special value of M 



B. Phase stability and phase diagram 



M S i 



7 

2 7 i 



(49) 



In order to outline the (t, r) phase diagram we shall 
use the information given above for the other two phases 
which have their own domains of stability in the (i, r) 
plane: N and FM. The FS stability conditions when 71 7^ 
become 



that is a solution of Eq. (43) for any value of r and 



t Si - --To 



11 
4 7l 2 ' 



(50) 



2 7 M-r-7iM 2 > 0, 



7M > 0, 



(45) 



(46) 



for which this procedure cannot be applied and should 
be considered separately. Note, that M$i is given by the 
respective horizontal dashed line in Fig. 4. The analysis 
shows that in the interval < t < 7 2 /2 the phase 
transition is again of first order; here 



3(1 - 7 2 )M 2 + 3771M + t - 7ir - 7 2 /2 > 0. (47) 



t {1) - - 



4(1 + 71) 2 ' 



(51) 



and we prefer to treat Eqs. (45) - (47) together with 
the existence condition fa 2 > 0, with </> given by Eq. (42), 
with the help of the picture shown in Fig. 4. 

The most direct approach to analyze the existence and 
stability of FS phase is to express r as of function of 
(M, t) from the equation of state (43), 



To find the equilibrium magnetization of first order phase 
transition, depicted by the thick line ACB in Fig. 4 we 
need the expression for equilibrium free energy of FS- 
phase. It is obtained from Eq. (12) by setting (w = 
0,v = 0) and substituting r, fa as given by Eqs. (41), 
(42) and (48). The result is 



J 



f (i) 
Jfs 



M 2 



2(M 7l - 7 /2) 2 



x <^ (l-7 2 )M 4 + 7 7l M 3 + 2 



«(1 - 7?) - 4r 



M z 



2 77l iM + i(i- -L.) 



(52) 



where M = M, 



eq ■ 



by the formula 



For the phase transition from N to FS phase (0 < t < 
7 2 /2), M eq is found by setting the FS free energy from the 
above expression equal to zero, as we have by convention 
that the free energy of the normal phase is zero. The 
value of Meq^ for positive t is obtained numerically and 
is illustrated by thick black curve AC in Fig. 4. When 
tg <t<0 the transition is between FM and FS phases 
and we obtain M$ from the equation f FS — f FM — 
(—t 2 /2), where /fm is the free energy of FM phase. The 
equilibrium magnetization in the above t-interval is given 



M (1) * 

AVJ -eq 



7 



2(1+ 7i) ^ 



(53) 



and is drawn by thick line CB in Fig. 4. 

The existence and stability analysis shows that for r > 
the equilibrium magnetization of the first order phase 



transition should satisfy the condition Jtfi 1 ' < M, 



M, 



(i) 



r (i) 

eq 



< 



By Mq 1 ^ we denote the positive solution of {M eq ) = 
and its i-dependence is drawn in Fig. 4 by the curve 
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with circles. Mm' is the smaller positive root of stability 
condition (47) and also gives the maximum of function 

r^(M); see Eq. (48). The function m£ } is depicted 
by the dotted curve AB in Fig. 4. When t S \ < t < 
t^g the existence and stability conditions are fulfilled if 



-t < M < Msi, where y—i is the magnetization of 
ferromagnetic phase and is drawn by a thin black line 
on the left of point B in Fig. (4). Here we have two 
possibilities: r > for y/^i < M < and r < for 

Af (1) < M < M S1 . To the left of t S i and t > t S2 , where 



We shall not discuss the region, t > 0, r < 0, because 
we have supposed from the very beginning that the tran- 
sition temperature for the ferromagnetic ordering T/ is 
higher then the superconducting transition temperature 
T s , as is for the known unconventional ferromagnetic su- 
perconductors. But this case may become of substantial 
interest when, as one may expect, materials with T f <T S 
may be discovered experimentally. 



Discussion 



tS2 



(54) 



the FS phase is stable and exists for Msi < M < y— t. 
Here r will be positive when Mg < M < and r < 
for M^ 1} > M > M S1 . When t < t S2 , M < ^f~t and r 
is always negative. 

On the basis of the existence and stability analysis we 
draw in Fig. 5 the (t, r)-phase diagram for concrete values 
of 7 and 71. As we have mentioned above the order of 
phase transitions is the same as for 71 = 0, see Fig. 3, 
Sec. III. The phase transition between the normal and 
FS phases is of first order and goes along the equilibrium 
line AC in the interval (t A = -y 2 /2 and t c = 0)- The 
function (t) is given by Eq. (48) with M$ from Fig. 
4. 

N, FM, and FS phases coexist at the triple point C 
with coordinates t = 0, and r^) = 7 2 /4(7i + 1). On 
the left of C for tg < t < the phase transition line 

of first order rig (t) is found by substituting in Eq. (48) 
the respective equilibrium magnetization, given by Eq. 
(53). In result we obtain 



*(*) 



7 



4(l + 7i) 



t. 



(55) 



This function is illustrated by the line BC in Fig. 5 that 
terminates at the tricritical point B with coordinates 
from Eq. (51), and 



(1) = 7 2 (2 + 7i) 
B 4(1+71)2' 



(56) 



To the left of the tricritical point B the second order 
phase transition curve is given by the relation (40). Here 
the magnetization is M = \f—i and the superconduct- 
ing order parameter is equal to zero (cf> — 0). This 
line intersects t-axis at t$2 and is well defined also for 
r < 0. The function ri X \t) has a maximum at the 
point (tsi, 7 2 /47i); here M = Msi- When this point 
is approached the second derivative of the free energy 
with respect to M tends to infinity. The result for the 
curves r^q (t) of equilibrium phase transitions (N-FS and 
FM-FS) can be used to define the respective equilibrium 
phase transition temperatures Tpg. 



The shape of the equilibrium phase transition lines 
corresponding to the phase transitions N-SC, N-FS, and 
FM-FS is similar to that of the more simple case 71 = 
and we shall not dwell on the variation of the size of the 
phase domains with the variations of the parameter 71 
from zero to values constrained by the condition 7 2 < 1. 
Our treatment from Sec. Ill of the magnetic susceptibil- 
ity tensor and the thermal quantities can be generalized 
in order to demonstrate the dependence of these quan- 
tities on 71. We shall not consider such problems. But 
an important qualitative difference between the equilib- 
rium phase transition lines shown in Figs. 3 and 5 can- 
not be omitted. The second order phase transition line 
r e (t), shown by the dotted line on the left of point B 
in Fig. 3, tends to large positive values of r for large 
negative values of t and remains in the second quadrant 
(t < 0, r > 0) of the plane (t, r) while the respective sec- 
ond order phase transition line ri (t) in Fig. 5 crosses 
the £-axis at the point t$2 and is located in the third 
quadrant (t < 0, r < 0) for all possible values t < tsz- 
This means that the ground state (at K) of systems 
with 71 = will be always the FS phase whereas two 
types of ground states, FM and FS, can exist for sys- 
tems with < 7 2 < 1. The latter seems more realistic 
when we compare theory and experiment, especially, in 
ferromagnetic compounds like UGe2, URhGe, and ZrZii2. 
Neglecting the 71 -term does not allow to describe the ex- 
perimentally observed presence of FM phase at very low 
temperatures and relatively low pressure P. 

The final aim of the phase diagram investigation is 
the outline of the (T, P) diagram. Important conclusions 
about the shape of the (T, P) diagram can be made from 
the form of the (t, r) diagram without an additional infor- 
mation about the values of the relevant material param- 
eters (a s , df,...) and their dependence on the pressure 
P. One should know also the characteristic tempera- 
ture T s , which has a lower value than the experimen- 
tally observed 0, S 0, B @ phase transition temperature 
(Tfs ~ IK) to the coexistence FS-phase. A supposi- 
tion about the dependence of the parameters a s and a f 
on the pressure P was made in Ref. HH. Our results 
for Tf 3> T s show that the phase transition temperature 
Tps varies with the variation of the system parameters 
(a s , a/, ...) from values which are higher than the charac- 
teristic temperature T s down to zero temperature. This 
is seen from Fig. 5. 
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FIG. 4: The dependence M(t) as an illustration of stability analysis for 7 = 1.2, 71 = 0.8 and w = 0. The parameters of the 
theory (r,t,-y, 71, w, . . . ) are defined by Eq. (11). The horizontal dashed lines represent the quantities M,$i given by Eq. (49) 
and Ms2 = 2Msi- The line of circles AS1S2 describes the positive solution of Eq. (48). The thick line AC gives the equilibrium 
magnetization for t > 0. The thick line BC represents the equilibrium magnetization for t < as given by Eq. (53). The dotted 
curve is the smaller positive solution of the stability condition (47). The thin solid line BS1S2 is the magnetization M = \/~^t. 
The arrow indicates the triple point C. A and B are tricritical points of phase transition. The point Si corresponds to the 
maximum of the curve (40) for t < 0, and the point S2 corresponds to r^\t) — in Eq. (40). 
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FIG. 5: The phase diagram in the (t,r) plane for 7 = 1.2, 71 = 0.8 and w — 0. The parameters of the theory (r, t, 7, 71, 
w, . . . ) are defined by Eq. (11). The domains of stability of the phases N, FM and FS are indicated. A and B are tricritical 
points of phase transitions separating the dashed lines (on the left of point B and on the right of point A) of second order 
phase transitions from the solid line ABC of first order phase transitions. The FS phase is stable in the whole domain of the 
(t, r) below the solid and dashed lines. The vertical dashed line coinciding with the r-axis above the triple point C indicates 
the N-FM phase transition of second order. 
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In systems where a pure superconducting phase is 
not observed for temperatures T ~ Tt or T ~ TfSj 
we can set T s ~ in Eq. (21). Neglecting T s in 
Eq. (21) and assuming that (T* /Tf) «lwe obtain that 
T c = Tps ~ {T*T f )V 2 . Note that the first (T*/T f )^ 2 - 
correction to this result has a negative sign which means 
that a suitable dependence of the characteristic temper- 
ature T* on the pressure P may be used in attempts 
to describe the experimental shape of the FM-FS phase 
transition line in the (T, P) diagrams of UGe2 and ZrZii2; 
see, for example ,Fig. 2 in Ref. Fig. 3 in Ref. Q, 
Fig. 4 in Ref. 5]. The experimental phase diagrams 
indicate that Tf(P) is a smooth monotonically decreas- 
ing function of the pressure P and Tf(P) tends to zero 
when the pressure P exceeds some critical value P c ~ 1 
GPa. Postulating the respective experimental shape of 
the function Tf (P) one may try to give a theoretical pre- 
diction for the shape of the curve Tps describing the FM- 
FS phase transition line. The lack of experimental data 
about important parameters of the theory forces us to 
make some suppositions about the behavior of the func- 
tion T*(P). The phase transition temperature Tp$ will 
qualitatively follow the shape of Tf (P) provided the de- 
pendence T*(P) is very smooth. This is in accord with 
the experimental shapes of these curves near the criti- 
cal pressure P c where both Tf and Tps are very small. 
The substantial difference between Tf and Tps a t lower 
pressure (P < P c ) can be explained with the negative 
sign of the correction term to the leading dependence 
Tps(P) ~ [T*{P)T f {P)] 1/2 mentioned above and a con- 
venient supposition for the form of the function T*(P). 

Eq. (21) presents a rather simplified theoretical result 
for Tc = Tps because the effect of M 2 \ip\ 2 coupling is not 
taken into account. But following the same ideas, used 
in our discussion of Eq. (21), a more reliable theoretical 
prediction of the shape of FM-FS phase transition line 
can be given on the basis of Eq. (40). Using the knowl- 
edge about the experimentally found shape of Tf (P) and 
the definition of the parameters r and t by Eq. (11) we 
substitute T = T FS (P) in Eq. (40). In doing this we 
have applied the following approximations, namely, that 
T s ~ for any pressure P, T FS (P C ) ~ T f (P c ) ~ and for 
substantially lower pressure (P < P c ), T f (P) > T FS (P). 
Then near the critical pressure P c , we easily obtain the 
transition temperature T F $ ~ 0, as should be. For sub- 
stantially lower values of the pressure there exists an 
experimental requirement (Tps — T s ) <C (Tf — Tps). 
Using the latter we establish the approximate formula 
(Tf — Tps) — 7 2 6y /l\OLf. The same formula for 
(Tf — Tps) can be obtained from the parameter ts2(T F s) 
given by Eq. (54). The pressure dependence of the pa- 
rameters included in this formula defines two qualita- 
tively different types of behavior of Tps (P) at relatively 
low pressures (P <C P c )- (a) T F s(P) ~ below some 
(second) critical value of the pressure (P' c < P c ), and 
(b) finite Tps(P) up to P ~ 0. Therefore, we can esti- 
mate the value of the pressure P' c < P c in UGe2 , where 
T F s(P' c ) ~ 0. It can be obtained from the equation 



T f (P' c ) = (l 2 b) /2 h 2 af) provided the pressure depen- 
dence of the respective material parameters is known. 
So, the above consideration is consistent with the theo- 
retical prediction that the dashed line in Fig. 5 crosses the 
axis r — and for this reason we have the opportunity 
to describe two ordered phases at low temperatures and 
broad variations of the pressure. Our theory allows also 
a description of the shape of the transition line Tps(P) 
in ZrZn2 and URhGe, where the transition temperature 
T F s is finite at ambient pressure. To avoid a misunder- 
standing, let us note that the diagram in Fig. 5 is quite 
general and the domain containing the point r = of the 
phase transition line for negative t may not be permitted 
in some ferromagnetic compounds. 

Up to now we have discussed experimental curves of 
second order phase transitions. Our analysis gives the 
opportunity to describe also first order phase transition 
lines. Our investigation of the free energy (12) leads to 
the prediction of triple (C)and tricritical points (A and 
B); see Figs. 3 and 5. We shall not dwell on the possible 
application of these results to the phase diagrams of real 
substances, where first order phase transitions and multi- 
critical phenomena occur; see, e.g., Refs. j^EI^, where 
first order phase transitions and tricritical points have 
been observed. The consideration of such problems, in 
particular, the explanation of the phase transition lines 
in Refs. |36ll39l | requires further theoretical studies, that 
can be done on the basis of a convenient extension of the 
free energy (12). For example, the investigation of vortex 
phases in Ref. [33 needs taking into account the gradient 
terms (4). Another generalization should be done in or- 
der to explain the observation of two FM phases [3(| |3^ • 
Note, that the experimentalists are not completely cer- 
tain whether the FS phase is a uniform or a vortex phase, 
and this is a crucial point for the orientation of the fur- 
ther investigations. But we find quite encouraging that 
our studies naturally lead to the prediction of the same 
variety of phase transition lines and multicritial points 
that has been observed in recent experiments [361 l39| . 



V. ANISOTROPY EFFECTS 

Our analysis demonstrates that when the anisotropy 
of Cooper pairs is taken in consideration, there will be 
no drastic changes in the shape the phase diagram for 
r > and the order of the respective phase transitions. 
Of course, there will be some changes in the size of the 
phase domains and the formulae for the thermodynamic 
quantities. It is readily seen from Figs. 6 and 7 that 
the temperature domain of first order phase transitions 
and the temperature domain of stability of FS above T s 
essentially vary with the variations of the anisotropy pa- 
rameter w. The parameter w will also insert changes 
in the values of the thermodynamic quantities like the 
magnetic susceptibility and the entropy and specific heat 
jumps at the phase transition points. 

Besides, and this seems to be the main anisotropy ef- 
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feet, the w- and v-terms in the free energy lead to a sta- 
bilization of the order along the main crystal directions 
which, in other words, means that the degeneration of the 
possible ground states (FM, SC, and FS) is considerably 
reduced. This means also a smaller number of marginally 
stable states. 

The dimensionless anisotropy parameter w = u s j (b s + 
u s ) can be either positive or negative depending on the 
sign of u s . Obviously when u s > 0, the parameter w 
will be positive too and will be in the interval < w< 1 
to ensure the positiveness of parameter b from Eq. (10). 
When w < 0, the latter condition is obeyed if the original 
parameters of free energy (3) satisfy the inequality — b s < 
u s < 0. 

We should mention here that a new phase of coexis- 
tence of superconductivity and ferromagnetism occurs as 
a solution of Eqs. (13). It is defined in the following way: 



1 



M 



1-71 L 
1 

1^1 



7l(*+^)-r 



(57) 



and 



2wsin (6» 2 - 0i ) = 7 M , cos (0 2 - 0i) ^ . (58) 

In the present approximation the phase (57)-(58) is un- 
stable, but this may be changed when crystal anisotropy 
is taken into account. 

We shall write the equations for order parameters M 
and <pj of FS phase in order to illustrate the changes when 



7l M 2 



(1-w) 



> 0, 



and 



(1-w - 7l 2 )M 3 
3 

±- lll M 2 - 



t(l 



w 



2 



(59) 

(60) 
0, 



where the meaning of the upper and lower sign is the 
same as explained just below Eq. (44). The difference 
in the stability conditions is more pronounced and gives 
new effects that will be explained further, 



(2 - w)jM - r - 7i M 2 



1 



w 



>0, 



7M -wr- wjiM 2 > 0, 



(61) 



(62) 



and 



3(1 - w - 7 2 )M 2 + 3771M + t(l -w)- 7 2 /2 - 7ir 
1 — w 



> 0. 
(63) 



The calculations of the phase diagram in (t, r) param- 
eter space are done in the same way as in case of w = 
and show that for w > there is no qualitative change 
of the phase diagram. Quantitatively, the region of first 
order phase transition widens both with respect to t and 
r as illustrated in Fig. 6. On the contrary, when w < 
the first order phase transition region becomes more nar- 
row but the condition (62) limits the stability of FS for 
r < 0. This is seen from Fig. 7 where FS is stable above 
the straight dotted line for r < and t < 0. So, purely 
superconducting (Meissner) phases occur also as ground 
states together with FS and FM phases. 



VI. CONCLUSION 

We investigated the M-trigger effect in unconventional 
ferromagnetic superconductors. This effect arises from 
the M^iV^-couphng term in the GL free energy and 
brings into existence a superconductivity in a domain 
of the system's phase diagram that is entirely occupied 
by the ferromagnetic phase. The coexistence of uncon- 
ventional superconductivity and ferromagnetic order is 
possible for temperatures above and below the critical 
temperature T s , that corresponds to the standard second- 
order phase transition from normal to Meissner phase - 
usual uniform superconductivity in a zero external mag- 
netic field which occurs outside the domain of existence of 
ferromagnetic order. Our investigation has been mainly 
intended to clarify the thermodynamic behavior at tem- 
peratures T s < T < Tf where the superconductivity can- 
not appear without the mechanism of M-triggcring. We 
have described the possible ordered phases (FM and FS) 
in this most interesting temperature interval. 

The Cooper pair and crystal anisotropies have also 
been investigated and their main effects on the thermo- 
dynamics of the triggered phase of coexistence is estab- 
lished. In discussions of concrete real material one should 
consider the respective crystal symmetry. But when the 
low symmetry and low order (in both M and ip) 7-term 
is present in the free energy, the dependence of essential 
thermodynamic properties on the type of crystal symme- 
try is not substantial. 

Below the superconducting critical temperature T s a 
variety of pure superconducting and mixed phases of co- 
existence of superconductivity and ferromagnetism exists 
and the thermodynamic behavior at these relatively low 
temperatures is more complex than in known cases of 
improper ferroelectrics. The case Tf < T s also needs a 
special investigation. 

Our results are referred to the possible uniform super- 
conducting and ferromagnetic states. Vortex and other 
nonuniform phases need a separate study. 

The relation of the present investigation to proper- 
ties of real ferromagnetic compounds, such as UGe2, 
URhGe, and ZrZn 2 , has been discussed throughout the 
text. In these compounds the ferromagnetic critical 
temperature is much larger than the superconducting 
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FIG. 6: Phase diagram in the (t, r) plane for 7 = 1.2, 71 = 0.8, and w = 0.4. The meaning of lines and points is the same as 
given in Fig. 5. 



20 




FIG. 7: Phase diagram in the (t, r) plane for 7 = 1.2, 71 = 0.8, and w = —2. The straight dotted line for r < indicates an 
instability of the FS phase. The meaning of other lines and notations is the same as given in Fig. 5. 
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critical temperature (Tf 3> T s ) and that is why the 
M-triggering of the spin-triplet superconductivity is very 
strong. Moreover, the 71-term is important to stabilize 
the FM order up to the absolute zero (0 K), as is in 
the known spin-triplet ferromagnetic superconductors. 
Ignoring |l5j | the symmetry conserving 71-term does not 
allow a proper description of the known real substances 
of this type. More experimental information about the 
values of the material parameters (a s ,a/,...) included 
in the free energy (12) is required in order to outline 
the thermodynamic behavior and the phase diagram 
in terms of thermodynamic parameters T and P. In 



particular, a reliable knowledge about the dependence of 
the parameters a s and a,f on the pressure P, the value 
of the characteristic temperature T s and the ratio a s /af 
at zero temperature are of primary interest. 
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